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We investigate the breaking of global statistical isotropy caused by a dark energy component with 
an energy-momentum tensor which has point symmetry, that could represent a cubic or hexagonal 
crystalline lattice. In such models Gaussian, adiabatic initial conditions created during inflation 
can lead to anisotropies in the cosmic microwave background whose spherical harmonic coeflicients 
are correlated, contrary to the standard assumption. We develop an adaptation of the line of sight 
integration method that can be applied to models where the background energy-momentum tensor is 
isotropic, but whose linearized perturbations are anisotropic. We then show how this can be applied 
to the cases of cubic and hexagonal symmetry. We compute quantities which show that such models 
are indistinguishable from isotropic models even in the most extreme parameter choices, in stark 
(N ■ contrast to models with anisotropic initial conditions based on inflation. The reason for this is that 

the dark energy based models contribute to the CMB anistropy via the inegrated Sachs- Wolfe effect, 
which is only relevent when the dark energy is dominant, that is, on the very largest scales. For 
inflationary models, however, the anisotropy is present on all scales. 
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^ ' I. INTRODUCTION 

Accurate measurements of the anisotropies in the cosmic microwave background (CMB) have lead to a revolution 
in cosmology. They have provided strong evidence for the basic tenets of the ACDM cosmology. This model has three 
^ \ matter/energy components: baryonic matter, cold dark matter and dark energy in the form of cosmological constant. 

. It also postulates a spectrum of near scale invariant, adiabatic perturbations, thought to be created during inflation, 
^-r ■ These are usually assumed to been to be compatible with Gaussianity and global statistical isotropy (GSI). 

. However, both these crucial properties have recently been questioned by detailed analyses of the data from the 
Wilkinson Microwave Anisotropy Probe (WMAP) [l|. Non-Gaussianity is expected at some level if an inflationary 
epoch is the origin of the adiabatic density perturbations. The present claims of a detection of non-Gaussianity based 
on the bispectrum are somewhat higher than predicted by standard slow roll inflation, but could be an indication 
I of non-standard physics during inflation. Violations of GSI, at first sight, appear to be much more worrying, and the 
analyses which have found such properties have often referred to them as "CMB anomalies" . 
. !^ " The first such property identified was the low CMB quadropole. This was first noted on the basis of observations 
, from the Cosmic Background Explorer (COBE) satellite 0, but its significance was only seriously questioned when 
' it persisted in the WMAP data. Recent results from the WMAP collaboration claims that it is still compatible with 
5^ \ the ACDM model Q- However, these initial claims have spawned a series of paper which have drawn attention to 
a number of unusual properties. These features include a near vanishing of the angular correlation function C(9) on 
scales 9 > 60° H, Q, asymmetries in the distribution of the power spectrum J?, d, Q, extreme cold spots [13, El and 
correlations between the spherical harmonic coeffficients |l2l. [isl. [l^. [isl. [l6l. [l7|. 

Obvious criticisms of these claims are: (1) there is some undentified instrumental or data processing effect which 
has been ignored in the analysis; (2) the subtraction of the galactic/extragalactic foregrounds has created artifacts in 
the maps; (3) the statistical significance has been overestimated. However, no such systematic effects have been found 
and many of the anomalies are seen at lower significance in the COBE data. Moreover they appear independent of 
frequency and the shear weight of the evidence seems to point to the fact that there is something unusual about the 
WMAP maps. 
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It seems, therefore, sensible to consider how the standard ACDM FRW paradigm with initial perturbations created 
during slow-roll inflation can be modified to provide an explanation for the observed anomalies. A range of ideas 
based anisotropic inflation have been proposed [H, [l^ which lead to an initial condition based explanation. A specific 
realization which has generated interest is due to Ackerman, Carroll and Wise ( ACW) [1^ , whose model results in a 
dipole anisotropy in the initial power spectrum which can be parameterized by 

P(k) = P(fc)(l+5.(k.n)2), (1) 

where k and n are unit vectors in Fourier and real space, and quantifies the level of anisotropy. The resulting CMB 
covariance matrix can be written in the form 

mi 7712 

1 + g*/i 

where the geometric coefficients 6imi£2m2 couple to £2 = {^1,^1 + 2} and mi to m2 — {mi, mi + f,mi + 2}. The 
can be computed from a modified version of a CMB code such as CMBFAST [2l| or CAMB The normal- 
ization factor in the denominator ensures that the cylindrical power spectrum ^FprSm^^™^™ ~ independent of 
the anisotropy parameter g^,. 

A recent analysis [2^ has used the WMAP V- and W-band maps to constrain the parameter g^,. They found 
statistically significant deviation from = with the best fitting value of g^, = 0.15 ± 0.04 and the preferred axis 
in the direction of {l,b) — (110°, 10°). We will use the apparent success of this model in explaining the data as a 
benchmark to compare with the alternative models which we will construct. However, we also note the recent work 



of ref. |24| which claims the ACW model (and more generally any inflationary model driven by a vector field) admits 
unstable solutions as perturbations cross the horizon. 

The oberved CMB fiuctuations are a convolution of the initial conditions and the transfer function which models 
the dynamical effects of the expansions of the Universe and the matter/energy components. An obvious alternative 
to an initial conditions origin for the observed anomalies is that they are due to anistropy in the transfer function. 
In particular, we explore the possibility that dark energy (with properties different to those of a pure cosmological 
constant) is responsible. Since dark energy has only made its presence felt since redshifts z w 1 (if the equation of 
state w w — 1, although there is also the possibility of some subdominant early dark energy if the equation of state 
is dynamical), it leaves an imprint on the CMB at the large angular scales required. This imprinting occurs via the 
Integrated Sachs Wolfe (ISW) effect, as photons experience time varying gravitational potentials along the line of 
sight to the surface of last scattering. 

The nature of dark energy is still a mystery (for a review of dark energy models, see ref. j25|). If dark energy is 
a regular perfect fluid, then the sound speed is equal to Cg = y/dP/dp = -^/w, where the equation of state w is the 
ratio of pressure P to energy density p. Since w < —1/3 to achieve the observed acceleration, this means that the 
sound speed would be imaginary, leading to instabilities in the fluid. This work is based on the postulate that dark 
ener gy is an elastic fluid, with non-zero anisotropic stress (models with anisotropic stress have also been investigated in 
refs. [26l[27l.[28t). The degree of elasticity is controlled by the shear modifli of the fluid. If these moduli are sufficiently 
large, the fluid can be stable even when the pressure is negative, making it suitable as a macroscopic model to describe 
dark energy. For example, a non-zero isotropic shear modulus p. modifies the sound speed to Cg = w + Ap/ [3(1 + 'w)p], 
and hence if p is large enough the sound speed is real. 

It is possible that a number of shear moduli could characterize the elastic properties of the fluid. In this case the 
dark energy will generate anisotropic perturbations, as the sound speed will be direction dependent (initial studies 
of dark energy with an anisotropic equation of state w have also been recently carried out [29l. Isoj). Our work here 
extends a proof of concept of this idea in ref. [Slj, in which we computed the evolution of cosmological perturbations. 
We showed that isotropic initial conditions could lead to anisotropy in the case of cubic symmetry. Here, we calculate 
the observational effects on large angle CMB fluctuations. 

The paper is organized as follows. In Section|TT]we outline the formalism, firstly reviewing the perturbation equations 
required to compute CMB anisotropics and showing how mode-mixing occurs between scalars, vectors and tensors. 
We then construct the covariance matrix, which is non-diagonal, as correlations exist between angular modes. We 
then outline an efficient method to compute this covariance matrix. In Section IIIII we present the numerical results 
of our computations, in the process comparing these with the ACW model. Finally, we provide a discussion of our 
results and some concluding remarks. 
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II. FORMALISM 



A. Stress-energy of an anisotropic medium with isotropic pressure 

The General Relativistic treatment of an elastic medium has been studied in detail by Carter and others [H, [13, S, 
[ssj . The primary goal of this work was to understand the dynamics of neutron stars. Recently, we have used a similar 
approach to study the dynamics of an isotropic elastic medium acting as dark energy in a flat FRW cosmology [11] . 
Here, we consider the anisotropic case in a flat FRW framework, following on from work in ref. [sH - The original 
manifestation of this model was studied in refs. [13, [H]- 

The properties of a general clastic medium are characterized by its energy density p, pressure tensor P'^'' and 
anisotropic stress, which is specified by the shear tensor J]^'^'"^. The first two of these quantities affect the properties 
of the background space-time and its perturbations, while the latter affects only the perturbations. The pressure 
tensor obeys the symmetry and orthogonality relations P'^'^ = P^^"^'> (where the brackets denote symmetrization with 
respect to the indices), P^'^u^, = 0, where u^^ = a~^(l,0, 0, 0) is the fluid flow vector, such that there are 6 free 
components of P^''. The shear tensor obeys the relations S^"'"" = S]('^'')('"^) = ^p'^m'-^ = 0, such that 20 

components specify the linear shear response to perturbations [39j. 

In this work we assume that the pressure tensor is isotropic, so that it is characterized by a single scalar P though 
pfj-u _ P'^P", with P = wp, where the projection tensor is 7'*'^ — g^^'^ + w^w'^ and the flow vector is normalized by 
u^Ufj_ = — 1. This condition restricts the free components of Yjf^^f"^ ^ and therefore the class of models which we study. 
The more general case of an anisotropic pressure tensor would need to be embedded in a Bianchi background, which 
is beyond the scope of this work. 

In the FRW synchronous gauge, the perturbed energy-momentum components of an elastic medium with isotropic 
pressure are given by [1^ 

5T\ = ip + P) {d.,C + h/2) , (3a) 
ST\ = -{p + P)e, (3b) 
ST', = -S',P {d,e + h/2) - S^/j + , (3c) 

where h is the trace of the metric perturbation and P = {p + P)dP/dp is the bulk modulus of the fluid. The vector 
is the spatial displacement of the fluid worldlines with respect to the background coordinates - the time part 
vanishes due to the orthogonality condition ^''u^j = 0. The evolution equation for is given by the conservation of 
energy-momentum V^T^i, = 0, 

{P + P) iC +'HC)- 3/3ne + d'ST', = , (4) 

where dots denote derivatives with respect to conformal time and Ti. is the conformal Hubble parameter. The time 
derivative is the velocity of the fluid. 

Here we discuss two classes of model with an isotropic pressure tensor but an anisotropic shear tensor. For our 
numerical results we focus on the first of these (elastic dark energy with a cubic shear tensor). We emphasize that 
these models are very different to the multiply-connected spaces which have recently been used in an attempt to 
explain some of the CMB anomalies (see for example ref. [i^ for recent constraints). Here, we have a flat background 
with unconnected topology, and an additional stress-energy component whose perturbations have point symmetry. 



1. Cubic symmetry 

A model with cubic symmetry automatically has an isotropic pressure tensor (4]| with a single degree of freedom, 
the pressure scalar, or equivalently the equation of state w. There are two degrees of freedom in the shear tensor. If 
we orient the cubic cell with respect to Cartesian co-ordinates, the non-zero components are 

4/iL/3 , (5a) 
-2pl/5, (5b) 
MT , (5c) 

where pi^ and arc the longitudinal and transverse shear moduli. When pi^ = /ix = M this reduces to the isotropic 
case. The particular realization of the space-filling cubic unit cell fixes the values of the shear moduli - the procedure 
for doing this is outlined in ref. [4]| . and involves computing the change in energy of the unit cell under spatial 
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FIG. 1: Unit cells of two shapes with cubic and hexagonal symmetry. 



transformations. For example, the unit cell of the simple cube shown in the left panel of Fig. [T]has /il/p = 1/6 and 
/ix/p = 1/3. These values apply when the energy density is proportional to the total surface area of the unit cell (a 
realization being a Nambu-Goto domain wall network with P = — (2/3) p), or when the energy density is proportional 
to the total edge length of the cell (for example a Nambu-Goto string network with P ~ —(1/3) p). Values for /iL and 
/iT are computed in ref. [4l| . 

In order to compute the sound speeds at which perturbations propagate, one computes the eigenvalues of the Fresnel 
tensor, Q'^'^ = H^f^'^VpVa + Pv^j^v^, where is the direction of propagation [4l|. For the cube, a wave moving parallel 
to the normal of the faces (for example, [1,0,0]) propagates at a speed ^ w + 4/iL/ [3(p + P)], and two other 
polarization modes propagate in orthogonal directions at speed = p.T/{p + P). In order to ensure stability and 
causality, such that < < 1, one finds a minimum requirement of p-l/ p > — 3w(l + w)/4, and an upper limit of 
Pt / p < I + w. Analysis of waves propagating in other directions gives identical limits with /iL Pt and pT — > p^. 
The values of p which we study will always respect — 3u;(l + u;)/4 = /I™'" < ph/p, Pt/p < /2'"^^ = 1 + w. 

In the cubic case the spatial components of the energy-momentum tensor are then 

ST' J - -5',il3 - 2pl/3) (9,e' + h/2) - 2ML(%f ^ + - 2S', , (6) 

where 5"*^ is the anisotropy source matrix, and is given by 

/ 9(,e^-f/i%/2 9(,e)+/i^/2\ 

S',=Ap\ d^ye^+hy,/2 9(,e)+/i^/2 . (7) 

V + /i^,/2 d^,e^+h%/2 J 

The difference in shear moduli Ap ~ /ix — Pl quantifies the amount of anisotropy. If Ap is non-zero the standard 
decomposition of perturbations into scalar, vector and tensor (SVT) modes is no longer valid, since S^j cannot be 
expressed in terms of a single scalar, vector or tensor quantity. This means that mode-mixing occurs and scalar initial 
conditions can excite vorticity and gravitational waves. 

The equation of motion for the fluid displacement vector ^' is given by 

(P + P)iC + nC) - 3l3ne - I3{d'dj^^ + d'h/2) - pLidWjC + 9'9jeV3 + d'h'j - d'h/3) = V' , (8) 
the component V'' sourcing anisotropy. This is related to 5*'^ via the Einstein equations, and is given by 

/ {dydy + d.d')^'' + d^{dy^y + d,^') + dyh^y + d'h^, \ 
V = 2 d^S'j = Ap {d^d^ + d^d^)£,y + dy{d^£,^ + 9^^^) -f- d'^hy^ -f- d'-hy^ . (9) 

\ (9^9^ + 9y9y)e + d^ld.^s;'' + a^c^) + -t- dyh% y 



2. Hexagonal symmetry 

Unlike a cubic system, the pressure tensor is not automatically isotropic in the hexagonal case. There are three 

degrees of freedom in the shear tensor. If we align the Cg symmetry axis with respect to the z co-ordinate, the 
non-zero components are 

S^^/, = E^/y = ^p/i , = 4/iA/3 , (10a) 

= -2p/i , = = -2p^/i , (10b) 

^%\=p, SV. =SV. = MB. (10c) 
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In the X — y plane the pressure and shear tensor are isotropic. This is a property of 2D hexagonal structures [4l| . 

In order to satisfy isotropy of the pressure tensor in 3D, the properties of the space-filling unit cell are not completely 
free. The unit cell of a simple hexagonal prism is shown in the right panel of Fig. [T] If the edges surrounding the top 
and bottom hexagonal faces have unit length, the length L must be fixed so the pressure tensor is isotropic. If the 
energy density is proportional to the total edge length (such that P = —(1/3) p), L must also be of unit length, and 
if the energy density is proportional to the total surface area (such that P = —(2/3) p) then L = 2\/3. In both cases 
one finds /i/p = 1/4, ha/p = 3/16 and hb/p = 1/3. 

As in the cubic case one can compute the eigenvalues of the Fresnel tensor to derive constraints on the values 
of p. A mode moving in the direction [1,0,0] propagates at v"^ — w + Ap/ [?>{p + P)]. with two other polarization 
states propagating in directions [0,1,0] and [0,0,1] at speeds v'^ = p/{p + P) and = pb/{p + P) respectively. 
Similarly, a wave in the [0, 0, 1] direction has a speed v'^ = w + Apx/ [3(p + P)] , with two polarization states with 
= pb/{p + P)- There is no fixed direction to derive lower limits on pB and upper limits on pA, so one should solve 
the eigenvalue solutions to ensures the solutions respect stability and causality. These conditions are satisfied for the 
isotropic hexagonal cell discussed above. 

In the hexagonal case the spatial components of the energy momentum tensor are given by ([6]), with /.jl —>■ p, and 
an anisotropic source term 



+ /i^ J2 

+ 

+ /i%/2 + /i%/2 

where A/xa = Pa — P and ApB — pb ~ p 



S'j - Aa^b I d\ye^ + hy,/2 I - ^ [5^, (9,e~ + /i.J2) + 5^,<5/ {dk^ + h/2)] , 

(11) 



B. Scalar- Vector-Tensor decomposition 

We now require a convenient basis to expand the perturbed Einstein and fluid equations. Scalar, vector and tensor 
(SVT) quantities can be expanded in terms of an orthogonal set of basis functions ef (a = 1 ... 3) in Fourier space. 
We set e] = ki, so that the fluid displacement vector is given by (see for example refs. [13, IHiHl for further details 
of this decomposition) 

e;=e(°)fc.+e^-i)e,? + e(+^)e?, (12) 

with analogous expressions for other vector quantities. Here (0) denotes the scalar component and (±1) the two 
vector modes. The metric perturbation in Fourier space is given by 

with ti^hY = k'^hfj = /i^ = 0. The vector and tensor components of the metric perturbation are then constructed 
from the basis vectors by 

hY = h^-'^ej + h^+'^el (14a) 
hi = (e?e^ - e?e^) + h^^'^ (e?ef + e^^) , (14b) 

with (±2) denoting the two tensor modes. For computational purposes we will use the particular choice of and e^ 
(which satisfy the orthogonality condition) given by 

I ( k \ if \ 

^ I. . \ - k . , hk ■ (15) 

V^? + fc^ V ; ^kl^kl V- {kl^kl) } 

The fluid equations of motion can now be obtained by taking the scalar product of ([8]) with the basis vectors. This 
is equivalent to decomposing SVT quantities in terms of eigenfunctions of the Laplacian as in refs. [12, [i^. However, 
in this decomposition SVT modes do not decouple due to the anisotropic source term. For the cubic case they are 
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given by (with analogous expressions for the hexagonal case) 

[p + + - 3/3Hp) + {13 + ifiL/S) (fc^e^o) + kh/2^ + ifiLkv 

where £, = The source vector for the cubic case is 

(^y ~l~ ^z)^^ kx{ky^^ ~\~ kz^^) kyhxy kzhxz 

-ikl + kl)£,y - ky{k,t + k,t) - k^hy, - k^hy, 

~{kl + kl)t - k,{k,t + ky^) - k^h 
In the synchronous gauge the Einstein equations are given by 
o2g% = -3H^ -Hh + did'h/2-didjh'^/2, 
2a^G\ = d,h-djh\, 
2a^G'o = djh'^-d'h, 



kyhzy 



2n-n^] s' 



hS' 



(16a) 
(16b) 
(16c) 

(17) 



(18a) 
(18b) 
(18c) 
(18d) 



Similarly, projecting the Einstein equations with the basis vectors gives the set of first order (constraint) and second 
order (evolution) equations in Fourier space: 

Constraint: 

nh/2-k^T] = inGa^Sp, 

kfi 47rGa2(p + P)P' , 

= 167rGa2(p + p)|(±i), 



(19a) 
(19b) 
(19c) 



Evolution: 



h + 2Hh-2k^ri = -2A'!TGa^SP, 
h + Qf] + 2H{h + Qi]) -2k^ri = -leTrGa^H^") 



where the sources are 



5p 
SP 

n(") 
n(-i) 
n(+i) 
n(-2) 
n(+2) 



dP 



dp 



5p, 



= 2pL [^ke"> + h/2 + Sjyj + 3ele{S} , 

= 2pL (fc|^-i) + + ^eleiS'j , 

= -2pLh'--^^ - 2 

= ~2pLh'^+^^ - 2 



=2 J 



3 J 



=2 J 



3 J 



(20a) 
(20b) 
(20c) 
(20d) 

(21a) 
(21b) 
(21c) 
(21d) 
(21c) 
(21f) 

(21g) 



We have modified the CAMB software [24I to include an anisotropic elastic component in an otherwise standard 
cosmology. This involves treating each of the (±) sources separately, and evolving SVT modes simultaneously due to 
the coupling (normally SVT modes decouple so can be computed independently). Furthermore, since the source term 
S^j is dependent on the direction as well as the magnitude of the k vector in Fourier space, we evolve equations at a 
discrete set of directions and magnitudes of k. We give more details of this procedure in Section III El 
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FIG. 2: Angular dependence in fe-space for the (0, — 1, +1, — 2, +2) cubic source functions (top-left, top-right, midddle-Ieft, 
middle right and bottom respectively), with A/i arbitrarily set to unity. We have used the equal area Aitoff projection. 



C. Angular dependence induced by anisotropy 

It is easy to see the angular dependence induced by the anisotropy by considering the projection of the source term 



S'^j with the basis vectors. For an initial scalar curvature fluctuation in the metric, htj = {kikj 
and hexagonal source terms have the form 
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Angular mode functions can then be defined by projecting the source term 



/ 2 2 



i z oini 
p3 3n cinit 
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e"e 



(22) 



(23a) 
(23b) 
(23c) 



These cubic source functions are displayed in an Aitoff projection in Fig. (5] It is interesting to note the morphology 
of the maps in each case. The (0, —I, —2) modes have the full set of cubic symmetries, the scalar term having obvious 
symmetry by inspecting the terms for S'^^^ (x (kxky)^ + [k^kzY + {^y^z)^- However, the mirror symmetry along the 
fc^-axis is broken for the (-1-1, -1-2) modes. No vector or tensor fluctuations are generated towards the centre of a cubic 
face or corner - these correspond to the direction of the maximal and minimal scalar fluctuations. 

The hexagonal source has two terms of different weight depending on A/.tA and A/xp. In all cases though, S^^^^ = 
non-zero sources arc shown in Fig. [3] and as expected, arc isotropic in the k^ — ky plane. 



8 




In order to see the effect of tlie anisotropic sources on tlie perturbation equations, we can write down a series 
expansion in conformal time r of the coupled differential equations. For concreteness, we consider a universe consisting 
of cold dark matter (c), photons (r) and an elastic fluid component (e) with w = —2/3 and a cubic shear tensor. We 
assume the photon moments with £ > 1 are zero, which is a good approximation until tight coupling breaks down. In 
the expansion terms below we suppress the CDM component for clarity (the hexagonal case has /iL ^ m)- 




where X = X/ p and uj^ = flxH^- One can see how anisotropy in the fluid induces anisotropy in the metric fluctuations, 
which in turn are a source of large- angle CMB temperature fluctuations via the ISW effect. Additionally, vector and 
tensor metric components are excited due to mode-mixing. 
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D. CMB covariance matrix 



The CMB temperature is given by r(x, 11,77) = "^iv) [1 + where the perturbation is expanded in terms 

of spherical harmonics Yf„i (fi), 



A(x, n, ?/) = ^ ^ a^m(x, 7])Yem{n) ■ 

Fourier transforming this expression gives the muhipole coefficients 

a£™(x,To) = ^ / d^^e^-^ j dri^ r,:„(n)A(k, n, r) . 
The hne of sight integral solution for the photon perturbations is 

/■To 

A(k,n,ro)= VCp(k)VF(")(k,n) / dre*-''r("^f)(k, r) , 

"'0 



(25) 



(26) 



(27) 



where x = k{TQ — t) and = k • n. The label p (between -2 and 2) represents the different varieties of initial 
conditions with corresponding to scalars, ±1 the two vector modes and ±2 the two tensor modes. The source 
function T("'^')(k, r) is that computed as a function of position in fc-space (not just k = |k|). It represents the 
response of the n-th mode to initial condition type p - all five modes will generally be excited by any particular initial 
condition, (^p is a random variable describing the properties of type-p initial conditions. 

The functions represent properties of the scalar, vector and tensor modes and are given by 



FW(k,n 
f(-i)(k,fi 
F(+i)(k,n 
i^(-2)(k,n 
i^(+2)(k,n 



= 1, 



„2„3 I „2„3\ 



Hence, we can deduce that 

afm(x. To) 

where 



The /£"'(a;,k) functions can be written in terms of linear combinations J^^!^ of the spin harmonics (as described in 
appendix E]) by 



(2^) 



4:)(x,k)= / df7fty,:„(n)e"''-"F(")(k,n). 



I d'ke'''-Cp(k)^'"drT("''')(k,r)4::)(x,k), 



(28a) 
(28b) 
(28c) 
(28d) 
(28e) 



(29) 



(30) 



t(") 



/t^(.,k)=W-H^j(«)(k). 



We will assume that the initial conditions are Gaussian random initial conditions and therefore 

(C;,(ki)Cp,(k2)) = (27r2)2pp^(fci)<5(3)(ki -k2)Jp,p, . 

One can now compute the correlation matrix 

Qimifsms = «mi(x,To)a£„„,(x,To)) = ^ J k^dkPp{k)A'^^^P^^^^{k), 



where 



\n1n2p 



(fc) = (_,)^.-l«.|/.-|n.| / d^2.j(-);(k)j(-) (k)A(:-P)(k)A(-'^)(k), 



(31) 
(32) 

(33) 
(34) 
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and 

A("-)(k) = r drT("-)(k, . (35) 

Jo [fc(To-r)]l"l 

In order to represent the covariance matrix, it will be convenient to replace {£, m) with a single index s = £{£+!) + 171, 
that is CsiS2 = Cfimi«2m2 for Si = £i{£i + 1) + mi and S2 = ^2(^2 + 1) + "^2- 

E. Approximate treatment 

The linearized calculation of section Hi CI suggests that for a small anisotropy coefficent (that is, small A/i's) we can 
write the transfer functions as 

A("'^)(k) = 5„p4'')(fc) +i?("''')(fc)5(")(%,0^), (36) 

where ~ i?^^^'^'' and B^, "^'^^ — B^i^^'^K This consists of an isotropic piece A'f\k), which is non-zero only for 

the initial condition type p, and an anisotropic piece B^^'^^ {k) which is non-zero for all modes. This approximation 
significantly reduces computational time, as only several evaluations of A^"'^^(k) are required at each position in k- 
space for each k = |k|. With knowledge of the source function S^'^\9f,,4'j^), this is sufficient to extract the amphtude 
of the isotropic and anisotropic components of the transfer function. From this, one can evaluate the integral p4p to 
compute the correlation matrix elements. 

To save further computational time, each time the covariance matrix is evaluated (for example, with a different set 
of cosmological parameters), the integrals over the anisotropic source functions will be the same. Therefore, one can 
pre-compute matrices of the form 

n^],"= = (_i)«i-i"ili^2-|n2| J dr!^J,^f,^^\*4m^^5("i)5("^) . (37) 

With some simple algebra the correlation matrix ()33|) can be expanded in terms of these geometric matrix elements, 
and the angular information p7p can be supplemented with amplitudes to give the contributions to the correlation 
matrix. For the case of cubic symmetry, these elements couple £1 to £2 = {£i,£i + 2,£i + 4} and mi to m2 — 
{mi, mi + 1, mi + 4, mi + 8}. The resulting matrix is sparse, and we make use of this for storage purposes. A naive 
approach would lead to the storage of ©(^max) elements, most of them zero, while the more efficient method only 
scales as C'(^max)- 

In doing this, one also gains some insight into the relative contribution between the isotropic and anisotropic 
components. In Fig. 2] we show the both the ISW and total ISW -I- last scattering surface (LSS) components of 
the scalar transfer function, for a model with w = —0.4, /It = (J-t/p = /i™'" and /11 = fJ-h/p = /J™'" -t- 0.1, where 
^min _ 0.18. The justification for using these values is discussed in the following section, but using |A/2| — 0.1 allows 
one to use the approximate treatment outlined here rather than the full numerical evolution. The covariance matrix 
to ^max ~ 40 can be computed in only several minutes with this method, compared to the full evolution, which takes 
^ 1 day. We find errors of less than 5% for any matrix element given the various w and | A/i| ~ 0.1 which we consider 
when compared to the full treatment that is significantly more time consuming. 

III. RESULTS 

We now compute the covariance matrix for anisotropic elastic dark energy (AEDE) models with a cubic shear 
tensor. Our baseline cosmological parameters are those listed for the best fitting ACDM model in ref. [i^, with 
elastic dark energy taking the place of a cosmological constant. We use a scalar comoving curvature perturbation as 
our initial conditions and we investigate models with w = —0.2, —0.4, —0.6 and —0.8. When we vary the value of 
w, we also adjust the value of the Hubble parameter, Hq = 100 /i kmsec^^ Mpc^^, to keep the CMB peak positions 
consistent with the ACDM model. We also keep the physical matter and baryon densities ^Ij^h^ and il\^h^ constant, 
which ensures the CMB peak height also remain fixed. In this way the only difference in the CMB correlation matrix 
between models occurs at low £. 

Dark energy effects the ISW contribution to CMB anisotropies in two ways. The first is by its effect on the expansion 
rate (and the induced decay of gravitational potentials as the total equation of state changes). The second is by the 
growth of imhomogenities in the dark energy itself [46l. [47|. in which the sound speed plays an important role. The 
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FIG. 4: Isotropic contribution to the scalar transfer function Ai(k) for the ISW component (left) and total ISW + last- 
scattering surface (LSS) components (middle). Also shown is the anisotropic component (right) which is generated only via 
the ISW effect. We use a model with w = -0.4, = M™'" and = m""" + 0.1, where ^i""" = 0.18. 

case of an isotropic clastic dark energy was discussed in ref. [s^. The situation here is similar, except we now have a 
anisotropic sound speed along with a source of large-scale vector and tensor perturbations. 

We find that the largest deviations from isotropy occur for the case w = —0.2. This is due to two reasons - as w 
approaches —1 dark energy perturbations are suppressed, and a higher value of w results in dark energy domination 
at earlier times. The anisotropic contribution to the transfer function then extends to smaller scales, and hence 
contributes to higher i in the covariance matrix. This value of w would be incompatible with recent SN observations 
(see for example ref. (48j). but serves as a baseline for the size of deviations from isotropy that can be generated in 
these models. We tried using even more extreme values oi w > —0.2, but in this case the size of the effect actually 
decreased. The reason for this is there is some overlap in the k range of the ISW and LSS components of the transfer 
function, which can be seen by examining Fig. HI This overlap results in partial cancellation of the isotropic-anisotropic 
cross term, which is strongest for cases with w > —0.2. 

For the shear moduli, the minimum and maximum values of /l which we considered were /I™'" = — 3u;(l -I- w)/4, 
^max — 1 due to stability and causality constraints. We obtained roughly the same results by interchanging 

/It and and found maximal deviation from isotropy by fixing cither one to its minimum value. For the other 
value, it was convenient to use |A/i| = 0.1 for several reasons: (1) This allowed the approximate treatment to be used, 
which was computationally much faster and sufficiently accurate; (2) \A/2\ = 0.1 corresponds to roughly the level of 
anisotropy found in domain wall lattices, which could be a possible realization of anisotropic dark energy (4l| : (3) the 
level of anisotropy does not increase significantly for |A/i| > 0.1. 

The reason for the last point can be understood by considering the effect of the dark energy sound speed on the 
scalar ISW source term. In the isotropic case, the difference between a model with Cg = and 0.1 is much larger than, 
say between Cg = 0.1 and 1 js^. In the anisotropic case, fixing one of the /x's to zero is equivalent to an anisotropic 
scalar sound speed which is zero in several spatial directions, and a maximal value of = 4|A/2|/[3(1 + w)] in other 
directions, with an appropriate interpolation between the two. The scalar source term is therefore rather insensitive to 
|A/i| > 0.1. A caveat here is that vector and tensor modes are sourced by the anisotropy, whose size are proportional 
to A/i. However, we find that the vector and tensor contributions to the covariance matrix decrease as a function of 
£ faster than the anisotropic scalar contribution, and the scalar terms dominate for all but the lowest £ modes. 

In the following discussion then, we will set /It = /I™™ and /2l = /I™'" -1-0.1. We will frequently compare the 
AEDE results with those of the ACW model. For this purpose, we will orient the preferred axis of the ACW dipole 
anisotropy along the equator, with 9^, = and 0* = 90° (see ref. [l^] for notation). In this case, the {m,m + 1} 
correlations of the covariance matrix vanish. 

The resulting numerical computation of the total CMB covariance matrix for £ = 2 ... 6 is shown in Fig. O for both 
the ACW and an AEDE model with w = —0.2. For visual clarity, we plot the expression 

{\Cs,sA - S,,.,Cf;Z) ^4^2 ih + 1) (^2 + l)/(27r) . (38) 

Notice that we use the norm of the covariance matrix, since this matrix can (and does in the AEDE case) have 
imaginary components, resulting from the i prefactor in eqn. (j34p . This is perfectly natural and this poses no problem 
for creating map realizations, since the matrix is Hermitian. We subtract the isotropic contribution of the covariance 
matrix in order to show the anisotropic components more clearly. The isotropic part can be computed by calculating 
the cylindrical Ci, defined by 



(39) 
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FIG. 5: Representation of the CMB correlation matrix (see text for details) in (j-iK)^ for the ISW component (top-left) and 
the total ISW + LSS (top-right), for an AEDE model with w = -0.2, fiT = /i™'" and /il = m™'" + 0.1. Also shown is the ACW 
correlation matrix for = 1 (bottom-left) and = 0.1 (bottom-right). 

These Ce's are shown for a range of models in Fig. [6l which we discuss in more detail below. The remaining £ factors 
in eqn. ([38|) are the analogous quantity to the usual £{£ + l)/(27r) when plotting the CMB power spectrum. 

Referring to Fig. [51 the thicker lines inside each plot correspond to the {2i+ 1) blocks for each m mode. Wc show the 
ISW and total ISW-I-LSS components for the AEDE model, and also compare these to the ACW covariancc matrix 
with = 1 and 0.1. There are several features to notice about these plots. Mode couplings occur in the ACW model 
between {£, i + 2} and in the AEDE model between {i, £ + 2} and {£, £ + 4}. These couplings are expected due to 
cubic symmetry. The planar angular dependence of anisotropy for the ACW model is also apparent, due to our choice 
of preferred axis. This results in more power in to = ±£ modes for each £, which is apparent along the main diagonal 
of the matrix. 

The size of the anisotropic terms from £ = 2 ... 6 arc comparable between the ISW only AEDE model with w = —0.4 
and ACW model with — 1. Furthermore, these terms are larger when considering the full ISW+LSS spectrum 
compared to the ACW model with = 0.1. However, one also notices a much faster decrease in the size of the terms 
as a function of £ compared to ACW, since anisotropy is only sourced at low £ in the AEDE model. 

Another important feature to notice is the change in amplitude of the anisotropic terms in the AEDE model when 
considering the full ISW+LSS spectrum. These are reduced by a factor of ^ 2 compared to the ISW only case, so 
one cannot simply add an isotropic Ce component - the cross term between the ISW+LSS is also important, as we 
discussed previously. The reduction in amplitude results from a partial cancellation in the isotropic-anisotropic cross 
term when compared to the ISW only case. 

Fig. [6] shows the cylindrical C^'s for a range of AEDE models with w — —0.2, —0.4, —0.6 and —0.8, again using 
p.T = M™'" and p.^ = yU™'" + 0.1 (except for w = —0.8 where we use Afi = 0.08 to comply with the causality constraint). 
This quantity is approximately the same for models with switched longitudinal and transverse shear moduli, which 
changes the sign of the anisotropy parameter Ap,. The effect of switching the sign of A/2 on the covariance matrix 
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FIG. 6: Cylindrical power spectrum Ce for w = —0.2 (top-left), w — —0.4 (top-right), w = —0.6 (bottom- left) and w — —0.8 
(bottom-right). The solid curve shows the ISW contribution for /it ~ /i™'" and fiL = ft™'" + 0.1 (the exception is for w = —0.8, 
where Ap. — 0.08 due to causality and stability bounds on the sound speed), and the dotted curve shows the total ISW+LSS 
component. Also shown for comparison is the cylindrical Cf for the ACW model (dashed curve). 

results in the {£,£ + 2} and {£,£ + 4} couplings remaining roughly constant, while the diagonal terms change sign 
(whilst keeping the Ce approximately constant). 

In anisotropic models it is known that the cosmic variance of the Cg is greater than an isotropic model with the 
same Ce [4§|. Physically, this is due a smaller number of degrees of freedom at each £ mode, that is the distribution is 
no longer X2e+i- T^^^ variance of the Cg can be found by diagonalizing the dimension (2£ + 1) square matrix for each 
£, which contains information on the m correlations. The result of this diagonalization, Chnem, can then be used to 
compute the variance using 



y-y f, 



(2^ + 1)2 '-'imhn 



(40) 



In Fig.[7]this is compared to the expected variance in an isotropic theory, o-^j^'^ — '2Cg/{2£+ 1). It is again apparent 
how the AEDE model approaches the isotropic limit much faster than the ACW model, but neither is particularly 
significant. 

Off-diagonal terms of the covariance (C^^, Cg^) are also non-zero in anisotropic models. We calculate this quantity 
numerically by simulating a large number of realizations 10^) of the sample covariance matrix. The result of these 
simulations for the AEDE and ACW models are shown in Fig. [H] For visual clarity, we plot the covariance relative to 
the isotropic value by 



(Cov(Q„ C,,) - Se,e, ^/fo) /^^^^^ ^2^iso 



(41) 



One can see that, relative the the main diagonal, the covariance (Ce-^jCi^) is larger for the AEDE model. This is 
apparent by examining the size of the £i — £2 correlations in Fig. [5l 

We now turn our attention to the question of whether the anisotropy is detectable by any observation. The optimal 
probability of distinguishing between two models with covariance matrices and , assuming model A is correct, 
in a cosmic variance limited experiment is given by 



In 



( P{{aen.}\A) \ 
{p{{ae,n}\B)J 



tr I 



(42) 
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FIG. 7: Fractional increase in cosmic variance over an isotropic model for (solid/dotted) the ISW/total contribution to the 
power spectrum, with w = —0.2, p,T = /i™'" and [ml = /i™'" + 0.1. Also shown is ACW result with = 1 (short-dashed) and 
=0.1 (long-dashed). 
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FIG. 8: (Left) Covariance between Ce for the ISW component of the power spectrum, using an AEDE model with w — —0.2, 
fiT = A™'" and fiL = m""" + 0.1. (Right) Same quantity for the ACW model with g^, = 1. 

This is equivalent to the expression used in ref. [Ill, but generalized to the anisotropic case; the derivation of this 
equation is given in appendix iBl In particular, we assume that the anisotropic model is correct, since we want to 
find the probability that it could be distinguished from the isotropic case. We also assume that both the isotropic 
and anisotropic model have the same C^, defined by eqn. (|39p . since wc wish to estimate the significance of only the 
anisotropic terms, not the overall power at each £. In this case, the first term on the right hand side is of eqn. (j42p is 
zero, so the probability is simply the ratio of the determinant of the covariance matrices. 

In Fig. [HI we show the both cumulative probability and contribution from each £ mode for AEDE and ACW models. 
The cumulative probability in the ACW model increases as roughly ^^j^x one would expect for an equal contribution 
from each mode, while the contribution from each i decreases markedly for i > 10 — 20 in the AEDE model. An 
extrapolation by eye suggests that a several a 'detection' would be possible for the ACW model with an ^max of several 
hundred, which is compatible with the results of ref. [1^. The AEDE models, however, do not reach this threshold 
when considering the aniostropy created by the total ISW + LSS components. 

Up to this point, we have only considered models with an scalar initial conditions. Certain models of inflation, 
however, are also expected to produce tensor fluctuations. The magnitude of these modes are model dependent, but 
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FIG. 9: Probability of distinguishing between an anisotropic and isotropic model each with the same Ce- The labeling is the 
same as for Fig. |B] where for the ACW model we use a value of ,g* = 0.1. On the left, we show the contribution from each i 
mode, and on the right the cumulative probability. 



power-law models of inflation, for example, have an observationally derived 2a limit on the tensor-to-scalar ratio of 
r < 0.42, at a pivot scale fco = 0.002 Mpc"^ [H]. 

The effect of isotropic elastic dark energy on tensor fluctuations was discussed in ref. [3^ . The additional source of 
anisotropic stress damps the evolution of tensor modes (which decay through horizon crossing), causing a reduction 
in CMB power on large-scales. In the anisotropic case, this damping will be dependent on the direction of the wave 
vector, causing correlations in the covariance matrix. Moreover, in the same way as vector and tensor modes are 
excited in the scalar case, scalar and vectors modes will be generated from purely tensor initial conditions. 

The probability of detecting anisotropy in AEDE models will change if tensor modes are present, so we have 
investigated this possibility using the same values of w as in the scalar case. Since the initial conditions for each of 
the scalar and tensor cases are uncorrelatcd, one can simply add the two covariance matrices. We find, however, that 
the detection probability does not change significantly for for several reasons. Firstly, the anisotropic contributions to 
the tensor covariance matrix along the main diagonal are actually opposite in sign to the scalar case - the nature of 
this can be seen by computing the individual components of eqn. (j37p . The tensor-tensor term, which is dominant in 
this case, distributes power along the main diagonal differently to the scalar-scalar term (for example, for the scalar- 
scalar term, more power goes to the m — modes). Secondly, the overall amplitude of the tensors are subdominant 
compared to scalars for observationally allowed values of r. Finally, we find that the probability of distinguishing 
between an anisotropic and isotropic model for the tensor covariance matrix also peaks around £ ~ 10. In summary, 
this fall off in anisotropic power on smaller scales appears to be a large stumbling block for this class models. 



IV. CONCLUSIONS 



We have developed a formalism for computing the CMB covariance matrix for models with an anisotropic dark 
energy component. It appears that the level of off-digonal terms which can be created by models where the pressure 
tensor is isotropic, but the perturbed energy-momentum tensor is not (via the elasticity tensor) is too small to be 
detected. This is true even in the most optimistic scenario where the values of w, /ix and /iL are chosen to give the 
largest possible signal. This is most explicitly illustrated by the fact that the dotted lines in Fig. [5] appear to top 
out a level below a 1-a detection when i^ax > 10. Note the contrast with the ACW model for which the curves 
increase (x ^^^x therefore inclusion of a sufficiently large number of ^-modes could lead to a detection as has been 
claimed [l^l . This is because the anisotropic impact of the dark energy only comes from the ISW effect which is only 
strong on the largest scales, whereas in the ACW model it is present on all scales. 
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This negative result is something of a disappointment since this idea was one of the few which could naturally lead 
to a non-diagonal covariance matrix from scalar density fluctuations created during inflation. Clearly for such a model 
to have a change of being detected, one needs to generate larger off-diagonal terms in the covariance model. This 
may be possible in more general dark energy models of this kind which have an anisotropic pressure tensor as well as 
an anisotropic elasticity tensor. This will require the model to be embedded in a more general anistropic spacetime 
such as Bianchi I. This possibility is presently under consideration. 
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APPENDIX A: CALCULATION OF ij"^ 
The standard expansion of a plane wave gives 

e-A' = ^z^(2^ + l)j,(x)P,(/.) = A7rJ2^'Mx)YUmLin) • (Al) 

I l,m 

This can be used in to show that 

4°) (x, k) = ATTi'j, (x)r,,„ (k) . (A2) 
We will use (0^, 4>j^) will refer to polar angles in fc-space. Using these angles we can write the basis vectors (fT5|) as 

ki = {sin 9 cos ipj^, sin 9 sin ipj^, cos j.) , 
ef = (sin0j,,-cos(/)^,O) , 



= (cos6^cos(pj^,cos9^sin(t)^, — sinOj.) . (A3) 
From these definitions, we can deduce that 

^ = , ^ = fc, sin + e\ cos , ( A4b) 

||=-fc., M = -e^os0,. (A4c) 
By computing derivatives with respect to 9 and </>, one can easily show that 

*™ IX sm 9^. d(l)f, X sm 9j, dcj)-^, 

= -^C - Ay,„(k) , (A5b) 

IX 39 f, X 39 



tJ' - i i i:4-'' - - ) (A5c) 



IX 



l^af - -^^ j r..(k) , (A5d) 

^ = -7^^- ^ = { -^4r) ^-(k) ■ (A5e) 

IX 89 1 x^ d9i, \sin9t d(hr. J 
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Using these expressions we can deduce that 



j(") _4_/-|n|M5) r(n),j^^ 



where the functions jl"'' are defined by 



r,™(k), 

1 a 



sinfl 

4m\^) = J7^^™(k), 



= sin^^j 



These can be written in terms of the spin-s harmonics 



Am ^(k) 



(^ + 1)! 



(£-1)! 
(£ + 1)! 



1/2 



(£-1)! 
(^ + 2)! 



1/2 



+iy£„(k)+-i Y^„^(k) 
2i 

-iy£,„(k) -_i Y«„,.(k) 



-,1/2 



(£-2)! 
(£ + 2)! 



+2l£rn(k) +-2 i^toi (k) 



^'^ ^r,™(k)-_2yf™(k) 



+2 J 



2i 



The spin-s harmonics can also be generated from the series expansion 



2£+ 1 (£ + to)! {l-m)\ 



-,1/2 



(A6) 

(A7a) 
(A7b) 

(A7c) 
(A7d) 
(A7e) 

(A8a) 
(A8b) 
(A8c) 
(A8d) 



47r (€ + s)! (€-s)! 
where {X^Y) is the binomial coefficient. 



„»m0 sin2^(6l/2) - s, r)(^ + s, r + s - m)(-l)^-'--^ cot2'^+-''-™(6'/2) , (A9) 



APPENDIX B: CALCULATION OF DISCRIMINATION EQUATION 

Given two models with covariance matrices C^-^^ and C*-^^, and the data vector a, the ratio of probabilities between 
models is 



\P{{a}\B)) 2 



^T.^(A) -1^ _ ^T.^(B) -1^ ^( ^£i^^ 

VdetC(^)y 



(Bl) 



We assume that C'"*' is the true covariance matrix so that a = AC^. where A{J(^y = is the Cholesky decomposition 
of the covariance matrix . For Gaussian initial conditions one has that {QCt') — ^w, so that 



In 



P{{a}\B) 



tr ( -^A] - tr ( AT'^C(^) -^A] + In 



det C^-^) 
detC(s) 



, /detC('4)\ 



(B2a) 
(B2b) 
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